
	if (ttime.gt.ttimeout) then
		nout=nout+1
	!	do j=1,8
	!	xg(1:3)=xg(1:3)-meshPoints(1:3,meshCellPoints(j,l))/8.0d0
	!	end do
!		if (mod(nout,25).eq.0) then
!			write(32) ttime,0,nCells,(Unew(1:3,i),pnew(i),Tnew(i),i=1,nCells)
!			if (rank.eq.0) write(*,'(a18,3e12.5)') 'ttime,dtmin,dt=',ttime,dtmin,dt
!		end if
	!	deltatprint=(ttimeout-(ttime-dt))
		do i=1,nCellsOut
		!	vtmp=(rhonew(meshCellOut(i))-rho(meshCellOut(i)) )/dt
			vtmp=rhonew(meshCellOut(i))
			Cellout(1:3,i)=Cellout(1:3,i)+Unew(1:3,meshCellOut(i))*(ttimeout-(ttime-dt))
			Cellout(4,i)=Cellout(4,i)+pnew(meshCellOut(i))*(ttimeout-(ttime-dt))
			Cellout(5,i)=Cellout(5,i)+vtmp*(ttimeout-(ttime-dt))
		end do
	!	write(*,*) rank,nCellsOut
		write(outFW) ttimeout,0,nCellsOut,((cellOut(j,i)/dtout,j=1,5),i=1,nCellsOut)
	!	if (rank.eq.0) write(*,'(a18,i10,3e12.5)') 'nout,ttimeout=',nout,ttimeout
		do i=1,nCellsOut
		!	vtmp=(rhonew(meshCellOut(i))-rho(meshCellOut(i)) )/dt
			vtmp=rhonew(meshCellOut(i))
			Cellout(1:3,i)=Unew(1:3,meshCellOut(i))*(ttime-ttimeout)
			Cellout(4,i)=pnew(meshCellOut(i))*(ttime-ttimeout)
			Cellout(5,i)=vtmp*(ttime-ttimeout)
		end do
		ttimeout=ttimeout+dtOut
	else
	!	if (rank.eq.0) write(*,'(a18,i10,3e12.5)') 'nout,ttimeout=',nout,ttimeout
		do i=1,nCellsOut
		!	vtmp=(rhonew(meshCellOut(i))-rho(meshCellOut(i)) )/dt
			vtmp=rhonew(meshCellOut(i))
			Cellout(1:3,i)=Cellout(1:3,i)+Unew(1:3,meshCellOut(i))*dt
			Cellout(4,i)=Cellout(4,i)+pnew(meshCellOut(i))*dt
			Cellout(5,i)=Cellout(5,i)+vtmp*dt
		end do
	end if
